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Abstract. Context. Hydrogen ionized by young, high-mass stars in starburst galaxies radiates radio recombination 
lines (RRLs), whose strength can be used as a diagnostic of the ionization rate, conditions and gas dynamics in 
the starburst region, without problems of dust obscuration. However, the lines are weak and only few extragalactic 
starburst systems have been detected. 

Aims. We aimed to increase the number of known starburst systems with detectable RRLs for detailed studies, 
and we used the line properties to study the gas properties and dynamics. 

Methods. We searched for the RRLs H91a and H92a with rest frequencies of 8.6 GHz and 8.3 GHz in the nearby 
southern Seyfert galaxy NGC 4945 using the Australia Telescope Compact Array with resolution of 3". This 
yielded a detection from which we derived conditions in the starburst regions. 

Results. We detected RRLs from the nucleus of NGC 4945 with a peak line strength integrated over the source 
of 17.8 mjy, making it the strongest extragalactic RRL emitter known at this frequency. The line and continuum 
emission from NGC 4945 can be matched by a model consisting of a collection of 10 to 300 HII regions with 
temperatures of 5 000K, densities of 10 3 cm -3 to 10 4 cm -3 and a total effective diameter of 2pc to 100 pc. The 
Lyman continuum production rate required to maintain the ionization is 6 x 10 52 s^ 1 to 3 x 10 53 s" 1 , which 
requires 2 000 to 10 000 05 stars to be produced in the starburst, inferring a star formation rate of 2Mgyr _1 to 
8Mq yr _1 . We resolved the rotation curve within the central 70 pc region and this is well described by a set of 
rotating rings that were coplanar and edge on. We found no reason to depart from a simple flat rotation curve. 
The rotation speed of 120 km s -1 within the central 1" (19 pc) radius infers an enclosed mass of 3 x 1O 7 M , and 
an average surface density with the central 19 pc of 25 000pc~ 2 , which exceeds the threshold gas surface density 
for star formation. 

Conclusion. We discovered RRLs from NGC 4945. It is the strongest known extragalactic RRL emitter and is 
suited to high-quality spectroscopic study. We resolved the dynamics of the ionized gas in the central 70 pc and 
derived conditions and star formation rates in the ionized gas. 
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1. Introduction 

NGC 4945 is a nearby almost edge-on spiral galaxy at 
a distance of (3.8 ± 0.3) Mpc ( |Karachentsev et al. (2007)] 
implying a scale of 19 pc per arcsec) . It has prominent dust 
lanes obscuring the nucleus, and is one of the brightest 
extragalactic sources seen by IRAS. The nuclear optical 
spectrum shows no sign of the Seyfert nucleus, and shows 



sorption (ATh ~2x 10 cm ) typical of Seyfert 2 galax- 
ies. The absorption was predicted to be much less at higher 



X-ray energies, and indeed at 100 keV Done et al. (1996) 



found it to be the second-brightest known Seyfert galaxy 
in the sky. X-ray imaging by Schurch et al. (2002) with 



a purely starburst emission (Krabbe et al. (2001) and ref- 



erences therein) . The presence of an active galactic nucleus 
(AGN) was finally settled when X-ray emission from a 
variable compact nuclear source was detected in the range 



X-ray Multi-Mirror-Newton and Chandra revealed a coni- 
cal plume extending 500 pc NW which they interpret as a 
mass-loaded superwind driven by the starburst, and they 
found a 6.4 keV iron Ka line and a Compton reflection 
component, which are characteristic of AGNs. 

A compact (< 30 mas) radio core with high bright- 



2keV to lOkeV (Iwasawaet al. (1993)) with heavy ab 



deceased 



ness temperature (10 7 K) detected by Sadler et al. (1995) 
at 2.3 GHz and 8.4 GHz indicates the presence of bright 
synchrotron emission, additional evidence for an AGN, 
though could possibly be radio supernovae. Mid-infrared 
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spectroscopy was carried out with the Infrared Space 



Observatory by Spoon et al. (2000) to penetrate the dust 
obscuration, and they found no high excitation lines like 
[Ne V] that are common in AGNs, implying huge obscura- 
tion of Ay > 160 assuming that a narrow-line region were 
present like that seen in, for example, the Circinus galaxy 
(which is similar in many ways to NGC 4945). 

The nucleus is the source of abundant infrared and 
molecular emission, and optical emission line images show 
gas outflowing from the nucleus perpendicular to the 
galaxy plane. This activity originates from a composite 
starburst and AGN. An obscured nuclear starburst ring 



with diameter 50 pc is seen in Paa (Marconi et al. (2000) ). 

The molecular content of the NGC 4945 nucleus 
is rich and varied, producing the strongest known ex- 
tragalactic molecular lines of many species. Over 80 
transitions from 19 molecules have been found with 
single dishes, principally OH, H 2 0, CO, HCN, HCO+, 
CH3OH, but also many others with the Swedish-ESO- 
Submillimetre Telescope (SEST) and Parkes from 
1.6 GHz to 354 GHz (|Gardner fc Whiteoak (1974J 



Whiteoak & Gardner (1974) 


Whiteoak et al. (1980) 


Whiteoak & Wilson (1990) 


Whiteoak & Gardner (1986) 



Dos Santos fc Lepine (1979)] |Batchelor et al. (1982) 



Henkel et al. (1990) 



Curran et al. (2001) 



Henkel et al. (1994) 



Wang et al. (2004)[ ). Many re- 
quire densities of 10 5 cm" 3 for excitation and yield cloud 
temperatures of T kin ~ 100 K. Wang et al. (2004) could 
determine isotope ratios of C, N, O, and S, and see 
isotope enrichment by ejecta from massive stars relative 
to that of fresh gas inflowing through a bar showing that 
the starburst is old enough to have affected the isotopic 
composition of the surrounding interstellar medium. 
Intcrferometric imaging of the HNC and HCO + molec- 



ular distributions by Cunningham & Whiteoak (2005) 
with the Australia Telescope Compact Array (ATCA) 
showed rotation of the molecular disk in the central 6" 
with rotation speed of 135 km s -1 . 

Imaging with very long baseline interferometry of the 



H2O masers by |Greenhill et al. (1997) resolved the rota- 
tion curve over the central 40 milliarcsec (0.7 pc diameter) 
with a velocity range of ±150 km s -1 inferring a central 
mass of 1.4 x 10 6 Mq. The rotation is in the same sense 
and in the same plane as the galaxy disk. 

CO has been imaged at resolution from 43" to 15" 
with SEST by |Dahlem et al. (1993)| |Ott et al. (200l)| 
and Mauersberger et al. (1996) and interferometrically 
at 4" with the Submillimeter Array (SMA) by 
Chou et al. (2007)| All groups found the CO to be 
strongly concentrated towards the centre in a disk 
of molecular material 16" x 11" (310 pc x 210pc). 



Chou et al. (2007) show that the disk rotates, within the 
central 5" (95 pc) radius in the plane of the galaxy disk 
with simple rigid-body circular rotation. At larger radius 
the isovelocity contours show an 'S' shaped asymmetry 
due to a bar potential, and at the centre there is an unre- 
solved kinematically decoupled component with a broad 
(340 km s -1 ) velocity range. They find that this last com- 



ponent is a good candidate for CO emission from the 
obscuring AGN torus, and too dense to be part of the 
starburst-driven molecular outflow. 

The ionized gas phase has been imaged in Ha 
+ Nil, [QUI], and Br 7 by |Moorwood fc Oliva (1994) 



Moorwood et al. (1996) and Lfpari et al. (1997) and in 



Paa by Marconi et al. (2000) The optical lines show a 



conical wind-blown cavity, and the infrared lines (Br 7 and 
Paa) show a starburst disk at the base of the cavity ex- 
tending over some 8" and so are embedded within the 
molecular disk which extends to 16". Despite the cone re- 
sembling Seyfert ionization cones, it lacks [O III] emission 
and so that origin has been excluded. Spectra in Ha + 
Nil show motions of ±500 km s _1 in the cone and opti- 
cal line ratios typical of low-ionization nuclear emission- 
line region galaxies. Infrared spectra covering Br7, Pf/3 



and shocked molecular hydrogen by Koornneef (1993) 
Moorwood fc Oliva (1994)| and |Spoon et al. (2003)| re- 
solve the nuclear rotating disk over the central ±250 pc, 
showing a velocity range of 500 km s _1 . The neutral hydro- 
gen optical emission region is within the molecular hydro- 
gen region, with shock excitation from a nuclear wind from 
the central cluster. The starburst seen in Br 7 is powerful 
enough to supply the 2 x 10 10 L Q radiated in the infrared. 

The ionized gas component can also be imaged using 
radio recombination lines (RRLs), which are not affected 
by extinction unlike at optical and near-infrared wave- 
lengths. In this paper we describe RRL imaging of the 
nucleus of NGC 4945 in which we detect H92a emission 
and resolve the rotation curve. RRLs have been used to 
derive mass, density and ionizing photon fluxes for the ion- 
ized gas in other galaxies, from which star formation rates 
have been derived. The diagnostic methods are described 
by Anantharamaiah et al. (2000)| and references therein, 
and in detail by Mohan (2002). The mere detection of 
RRL emission requires the presence of thousands of H II 
regions or of stimulated emission since a single H II region 
like the Orion nebula placed at that distance would be 
undetectably weak. 

There are now 15 known extragalactic RRL detec- 
tions (including NGC 4945), all in bright starburst galax- 
ies, which are listed by Roy et al. (2008) Most of those 



detections resulted from improving the search sensitivity 
by a factor of ten during the 1990s and by making surveys 
of promising bright candidates using the Very Large Array 
(VLA) and ATCA. Our observation of NGC 4945 reported 
here was the third detection made in our ATCA survey for 
H92a emission. The other two detections from this sur- 
vey (NGC 3256 and the Circinus galaxy) were reported by 



Roy et al. (2005)] and |Roy et al. (2008)] ). NGC 4945 has 
proven to be the strongest known extragalactic RRL emit- 
ter on the sky. 

We give velocities in the heliocentric frame using the 
optical velocity definition throughout this paper. 
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2. Observations 

We observed NGC 4945 with the ATCA simultaneously 
in the lines H91a and H92a with two orthogonal linear 
polarizations. The observing parameters and results are 
summarized in Table 1. 

Calibration and imaging were done using the AIPS 
software using standard methods. The flux-density scale 
assumed that PKS B1934-638 had a flux density of 2.99 Jy 
at 8295 MHz and 2.87 Jy at 8570 MHz, based on the 



-49 27 



Baars et al. (1977) flux-density scale. A phase calibra- 



tor was observed every half hour to correct the instru- 
mental phase response. A bandpass calibrator was ob- 
served every few hours for correcting the instrumental fre- 
quency response (bandpass). Phase corrections obtained 
from self calibration of the continuum source were ap- 
plied to the spectral line data. Continuum emission was 
subtracted from the line data using the method UVLSF 



(Cornwell, Uson & Haddad (1992) I in which the contin- 
uum is determined for each baseline by a linear fit to the 
spectrum. The final continuum and line images were made 
using natural weighting of the («, v) data to achieve max- 
imum possible signal-to-noise ratio. 

Uncertainties on the flux densities have an 11% rms 
random multiplicative component due to flux-density 
bootstrapping and atmospheric opacity, a 0.21 mJy rms 
random additive component due to thermal noise in a 
1 MHz channel and a systematic multiplicative component 
of 11 % rms due mainly to the uncertainty in the Baars et 
al. flux-density scale. 

3. Results 

The ATCA continuum and line images, and integrated 
spectrum are shown in Fig. 1, and the velocity field and 
position-velocity diagram are shown in Fig. 2. 

The continuum image shows a well resolved struc- 
ture extended along the plane of the galaxy in posi- 
tion angle 47°. Line emission was detected from the nu- 
clear region with a well resolved structure that is ex- 
tended along the plane of the galaxy like the contin- 
uum, though more compact and clumpy than the con- 
tinuum emission. The line emitting region has a total 
deconvolved extent of 8.1" x 1.7" (150 pc x 30 pc). The 
peak of line emission is coincident with the peak of con- 
tinuum emission within 1". The total area of line emis- 
sion is 1.5 times the naturally-weighted beam area (8.6 
times the uniformly- weighted beam area), or 4500 pc 2 . 
The H91a + H92ev spectrum integrated over the line 
emitting region shows a strong line detection with sig- 
nificance of 52 a, with centroid at 581kms _1 in good 
agreement with the CO systemic velocity of 585 km s" 1 



(Chou et al. (2007)1 and with the velocity of the HI ab- 
sorption towards the nucleus of 585 km s , though signif- 
icantly higher than the HI systemic velocity inferred from 
the whole galaxy of 557kms" 1 ( |Ott et al. (2001)[ ). The 
line FWHM is 220 km s -1 after deconvolving the instru- 
mental velocity resolution of 42kms~ 1 . 
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Fig. 1. Top: ATCA 8.3 GHz + 8.6 GHz continuum im- 
age with uniform weight of NGC 4945. Middle: ATCA 
uniformly weighted zeroth-moment image of H91a + 
H92q emission after continuum subtraction. Beamsize is 
1.4" x 1.2". Bottom: ATCA H91a + H92a line profile inte- 
grated over the line-emitting region in the zeroth-moment 
image. RMS noise is 0.16mJy beam -1 channel" 1 and the 
channel width is 36kms _1 . The region of integration is a 
box of size 10.0" x 10.4" centred on RA 13 05 27.48 dec 
-49 28 04.0 
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Table 1. ATCA Observational Parameters and Results for NGC 4945. 



NGC 4945 



Array configuration 
Date of observation 
Integration time 
Transitions 
l/ rest of RRL 

^band-centre 
Kyatemic , hclio , optical 

Distance 

Bandwidth, channels, IFs 

Spectral resolution 

Minimum baseline 

Beam size (natural weight) 

Beam size (uniform weight) 

Phase calibrator (B1950) 

Bandpass calibrator (B1950) 

Line and Continuum Properties 

Peak continuum flux density 

Total continuum flux density 

Noise per image channel 

Peak line flux density 

Integrated line flux 

Measured line width (FWHM) 

Deconvolved line width (FWHM) 

Ccntroid hclio. optical vcl. 

No. of beam areas where line is observed 



750D (1993jul25), 1.5D (1994scp21), 6C (94oct03) 
1993jul25, 1994sep21, 94oct03 
22. 7h 

H92q & H91a 

8309.38 MHz & 8584.82 MHz 

8295.00 MH z fe 8570.00 MHz 

560kms _1 jPaglione (1997)) 

3.8 Mpc 

64 MHz, 64, 4 

36 km s~ 1 

0.86 kA 

3.4" x 2.7" at -55° 
1.35" x 1.18" at -4° 
1320-446 

1226+023, 1320-446, 1934-638, 2251 + 158 



(538 ± 80) mjybeam 
(1424 ± 220) mjy 
0.21 mjy beam -1 
(17.8 ± 2.8) mjy 
(11.3 ± 1.8) X 10~ 22 Wm 
(225 ± 15) kms -1 
(221 ± 15) kms -1 
(581 ± 15) kms -1 
1.5 (naturally weighted) 



(natural weighting) 



The position-velocity diagram (Fig. 2) shows a well 
resolved symmetric structure indicating undisturbed ro- 
tation with a flat rotation curve extending in close to the 
nucleus. The rotation curve is modelled in Sec 5. 

4. Modelling Line Formation in the Ionized Gas 

Two types of models have been considered for the RRL 
emission from the nuclei of external galaxies: one based 
on a uniform slab of ionized gas and the other based on a 
collection of compact H II regions. Such models have been 



discussed by Anantharamaiah et al. (2000) and references 



therein, and are documented in detail by Mohan (2002) 



These models take as constraints the integrated RRL 
strength at one or more frequencies, the radio continuum 
spectrum, and the geometry of the line emitting region. 

For modelling the RRL emission from NGC 4945 we 
used the integrated line strength (17.8 mjy), line width 
(280 km s _1 ), size of the line-emitting region (equivalent 
to a 70 pc diameter sphere), continuum emission (1424 
mjy), and spectral index (-0.75) to constrain conditions in 
the ionized gas. Using the collection of HII regions model, 
models with 10 to 300 H II regions, all with T c - 5000 K, 
n c ~ 10 3 cm" 3 to 10 4 cm" 3 and a total effective diameter 
of the line-emitting gas of 2pc to 100 pc produced good 
matches to the line and continuum emission. Parameters 
derived for typical allowed models are given in Table 2. 

The range of possible filling factors can be estimated by 
comparing the total volume of the line-emitting gas (8 pc 3 
to 1 x 10 6 pc 3 derived from the total effective diameter) 
to the volume of the line-emitting region (5.3 x 10 5 pc 3 , 
assuming cylindrical geometry with diameter of 150 pc and 
height of 30 pc). The range of possible values is then 2 x 
10" 5 to 1. 



The inferred mass of ionized gas is 2 x 10 5 Mq to 
6 x 10 5 Mq, depending on the model conditions, which 
requires a Lyman continuum flux of 6 x 10 52 s _1 to 3 x 
10 53 s _1 to maintain the ionization. This flux is equivalent 
to the Lyman continuum output of 2000 to 10000 stars of 
type 05, which infers a star-formation rate of 2M yr _1 
to 8 Mq yr _1 when averaged over the lifetime of OB stars. 

This can be compared to star formation rates derived 
from other indicators following Hopkins et al. (2003) . 
Taking the peak 1.4 GHz flux density of 4.0Jybeam _i 
in the 30" x 18" beam of the ATCA at the nucleus of 



NGC 4945 (Elmouttie et al. (1997)1 yields a 1.4 GHz lu- 



minosity of 7.0 x 10 Wm and a corresponding star 
formation rate of 3.9 Mq yr^ 1 . This can be attributed to 
the central 1.8" (30 pc) diameter region following the 



argument of Elmouttie et al. (1997)| based on the high 



resolution 8.4 GHz image. The IRAS 60 /im and 100 /im 
flux densities yield a far-infrared (FIR) star formation 
rate of 4.1M yr _1 . These estimates agree well with the 
star formation rate estimated from the RRL emission of 
2M Q yr~ 1 to 8M Q yr _1 . The Ha and U-band-based star 
formation rates were not estimated due to the extreme 
absorption in the optical band. The supernova rate esti- 



mated by |Forbes fc Norris (1998)| from the 4.8 GHz flux 
density is 0.23 yr _i , and from the [Fe II] flux is much 
lower at 0.005 yr" 1 . These can be converted into star- 



formation rates following Fukugita & Kawasaki (2003) 
yielding rates of 19 Mq yr^ and 0.4MQyr -i . These 
bracket the range we estimated from the RRL emission. 

Our RRL-based SFR estimate assumed no dust ab- 
sorption between the OB stars and the HII regions 
producing the RRL emission. The effect of dust can 
be estimated by comparing the bolometric luminosity 



A. L. Roy et al.: Recombination Lines from NGC 4945 
Table 2. Derived properties for NGC 4945, using model results for T c = 5 000K 



5 



NGC 4945 



Electron temperature 
Electron density 
Effective size a 
Total ionized gas mass 

JVLyC 

No. 05 stars 



5 000 K 

1000 cm -3 to 10 000 cm -3 

2pc to 100 pc 

(2 to 6) xlO 5 M Q 

6 x 10 52 s _1 to 3 x 10 53 s" 
2 000 to 10 000 



The total size, ie (total volume of all H II regions) 



1/3 



(X D o1-rrl) from the stellar population needed to ionize 
the RRL-emitting gas to the observed FIR output (Lfir) 
from the region. If ionizing photons are absorbed by dust 
before causing ionization, then Lboi-RRL will be less than 
Lfir- The ratio Lboi— rrl / £fir thus gives the ratio of 
the number of photons absorbed by the gas to the number 
of photons absorbed by the dust. We derived £t>oi— rrl by 
taking our RRL-based I/LyC and dividing by the ^LyC / 
Lbol ratio of 0.29 calculated for a stellar population with 
Salpeter IMF, lower- mass cutoff of 1 Mq and upper-mass 
cutoff of 87 Mq, and using the mass- luminosity relation for 



OB stars in Table 5 of Vacca et al. (1996) This resulted 



in iboi-RRL = 1.2 x 10 y L Q to 5.7 x 10 y L Q , compared to 
Lfir from region of 1.6 x 10 10 L o (Brock et al. 1988). The 
ratio iboi-RRL / £fir shows that the gas volume in the 
nuclear region is illuminated by ionizing radiation that 
represents only 7.5 % to 35 % of the ionizing flux from the 
stars needed to power the FIR output; the rest, we assume, 
is absorbed by dust within the HII regions. Converting to 



opacity using Fig 2 of Petrosian (1972) yields dust opaci- 
ties of 1.3 neper to 3.4 neper within the HII regions. 

Given so much absorption, it is unexpected that the 
star-formation rate inferred from the RRL emission, which 
assumed zero absorption, should agree so well with that 
inferred from Lfir- An important assumption made while 
deriving the RRL-based star formation rate from the 
present inferred number of OB stars was that the stars 
are being formed at a steady rate that would maintain 
the present numbers. 



5. Dynamical Modelling of the Ionized Gas 

We used the rotation curve to constrain the gas kinemat- 
ics by fitting to the data a simple model consisting of a 
set of rings, coplanar, edge-on, with an initially flat ro- 
tation curve. The brightness of each ring was determined 
by deprojecting the observed zeroth-moment image to de- 
rive the radial distribution of the line intensity. The radial 
distribution showed a central peak and a ring of emission 
2.5" (50 pc) from the nucleus. 

We refined the model iteratively, varying the systemic 
velocity, rotation velocity, and velocity dispersion until we 
achieved a reasonably close match to the data, and found 
no reason to depart from a simple flat rotation curve. The 
biggest residual between the model and data was at the 



nucleus and a smaller one at the 50 pc ring. To improve 
the model, we refined the radial profile of line emission 
strength at the nucleus and the ring. The data required 
us to boost the central peak to 25 times the brightness of 
the more extended emission and to place that gas at the 
systemic velocity. This could correspond to gas on circular 
orbits moving transverse to the line of sight in front of 
the nucleus. The data also required us to put half that 
strength in the innermost ring (0.4" = 7.6 pc) with our 
best-fit rotation velocity. For the 50 pc ring, we tweaked 
the radius and brightness. The result was a very good 
match to the observed data, with residuals at 1.5 times 
the rms noise in the position- velocity diagram. 

The final model had a flat rotation curve with 
^systemic = 581 km s~\ ^rotation = 120 km s -1 , and 
^dispersion = 15 km s . We did not need to invoke a bar 
or radial motion, though the data do not strongly exclude 
cither. 

The rotation velocity of 120 km s _1 within the central 
1" (19 pc) radius infers an enclosed mass of 3 x 10 7 Mq. 
The water masers (Greenhillet al. (1997)1 infer an en- 
closed mass of 1 x IO^Mq within 0.3 pc radius, and so 
most of the 3 x 10 7 M Q is extended between 0.3 pc and 
19 pc radius of the nucleus. The average surface density 
within the central 19 pc is 25 000 Mq pc~ 2 , which exceeds 
the threshold gas surface density for star-formation of (3 
to 10) M Q pc~ 2 (Kennicutt (1989) I by four orders of mag- 
nitude. 

The 50 pc ring might be a Lindblad resonance, 
which would infer a bar and then our assump- 
tion of pure rotation would become invalid. The 
ring and nucleus may be seen also in CO spectra 



( |Whiteoak k Wilson (1990)[ |Dahlem et al. (1993)] ), and 



in OH absorption ( |Whiteoak k Wilson (1990)[ ). 

The RRL emission contains a bright peak on the nu- 
cleus at the systemic velocity. The continuum emission is 
smooth and does not show a similar peak and so we spec- 
ulate that this might be stimulated emission along the line 
of sight for which the gas velocity is transverse to the line 
of sight and so has a long maser gain path. 

6. Comparison with literature 



H I kinematics were studied by |Ott et al. (2001) with the 
ATCA at 3.5" resolution. They saw HI in emission over 
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the disk, with an extent of 22' x 4' (43kpc x 7.8 kpc) 
and in absorption towards the nuclear continuum emis- 
sion. They see an asymmetry in the large-scale H I rotation 
field, which they argue is due to the bar potential. The HI 
absorption towards the nucleus shows a uniform velocity 
gradient across the nucleus indicating solid-body rotation 
peaking at ±100 km s -1 either side of the systemic veloc- 
ity out to ±5" either side of the centre, with a sign of 
flattening beyond ±5" to ±7". For comparison, our H91a 
+ H92a rotation curve remained flat at ±120 km s" 1 ei- 
ther side of systemic all the way into the central 1". This 
difference could be due to the two-times lower spatial reso- 
lution of the HI observation compared to that of the H91a 
+ H92a observation (3.5" compared to 1.4"). Were we to 
smooth our position-velocity diagram, or had the noise 
in our position-velocity diagram been a little higher, the 
small wings that indicate a flat rotation curve would be 
less visible, giving the impression of solid-body rotation. 
Also, in edge-on systems with a flat rotation curve, one 
finds the density distribution often gives the appearance 
of solid-body rotation. 

Rotation of the larger-scale disk (> ±100") 
derived from CO(1-0) observations with SEST by 
Dahlem et al. (1993) at 43" resolution shows evidence for 
gas inflow that is consistent with the inflow model pro- 



posed by Abies et al. (1987) from HI. This was seen also 
in CO(2-iyby[Ott et al. (2001)] with SEST at 23" resolu- 



tion. However, inflow was not seen in the H91a + H92a 
dynamics within 4" of the nucleus. 



Isovelocity contours of CO by Chou et al. (2007) show 



an S-shaped asymmetry due to a bar potential at distances 
larger than 5" from the nucleus, and show simple circular 
rotation within that radius. The isovelocity contours of the 
H91a + H92a radio recombination line from the ionized 
gas component seen within 4" radius of the nucleus like- 
wise shows simple circular rotation with no asymmetry, 
consistent with the CO on that scale. 

The rotation of the nuclear disk derived from CO(3- 
2) observations with SEST by |Mauersberger et al. (1996) 



with 15" resolution and with SMA by |Chou et al. (2007) 
at ~ 4" resolution show that within ±5" of the nucleus, 
the velocity changes from 390 kms" 1 to 730 kms -1 , which 
is broader than the velocity range of the H91a + H92a 
emission (460 km s" 1 to 700 km s" 1 ) over the same region. 
The CO emission seen by both 

and 



Mauersberger et al. (1996) 



Chou et al. (2007) 



shows two local intensity peaks, with a radius of 6" to 8" 
from the nucleus. The H91a + H92a line emission radial 
distribution showed a central peak and a ring of emission 
with radius 2.5" (50 pc) from the nucleus. Thus, the ring 
in ionized gas lies within the ring in the molecular gas, 
consistent with the illustration by Spoon et al. (200*3)] 
Fig. 8. 

The dynamical mass within the central 100 pc ra- 



dius derived from the CO by Mauersberger et al. (1996) 



is 8 x 10 M© and the dynamical mass within 0.3 pc radius 



derived from the water masers by Greenhill et al. (1997) 



central 1" (19 pc) radius derived from the H91a + H92a 
line is 3 x 1O 7 M0, which lies between the masses within 
the larger and smaller radii, as expected. 

Rotation of the nuclear ionized gas is traced also by 



Pf/3 by Spoon et al. (2003) They find ionized gas within 
the central ±5" around the nucleus from the same re- 
gion in which we find H91a + H92a emission. Comparing 
their position-velocity diagram of Pf/3 (their Fig. 7) to 
our position-velocity diagram of H91a + H92a (our Fig. 
2), the diagrams look essentially the same, having the 
same central velocity gradient, the same rotation rate, the 
same systemic velocity, and the same overall extent. Thus, 
both emission lines are tracing the same gas, as expected. 



Spoon et al. (2003) propose a geometry in which the HI 



Pf/3 originates from the central ionized and rotating disk 
(their Fig. 8). The H91a + H92a kinematics supports this 
picture. Similar rotation curves are seen in HCN, HNC, 



and HCO+ by Cunningham fc Whiteoak (2005) though 



the molecular material is seen more strongly concentrated 
into an edge-on ring with weak emission or absorption to- 
wards the nucleus, unlike the ionized gas which is always 
in emission. 



Chou et al. (2007) identify a kinematically decoupled 



component in their CO position-velocity diagrams along 
the minor axis, showing broad velocity width at or near 
the disk centre, as a strip in velocity at the origin 
of the position- velocity diagram spanning 400 kms" 1 to 
750 km s" 1 in their Fig. 5. They identify this g as as a 
promising candidate for the circumnuclear molecular torus 
invoked by AGN unification models. For comparison, we 
show the position-velocity diagram in H91a + H92a along 
the minor axis in Fig. 3. The H91a + H92a shows a broad 
velocity range at the position of the nucleus (panel 252), 
spanning 430 kms -1 to 750 kms" 1 , similar to that seen 



in CO by Chou et al. (2007) This broad velocity width 
in the centre of our position-velocity diagram along the 
minor axis underlines that the rotation curve is flat; at 
this position close to the centre, one sees the full reced- 
ing and the full approaching velocities, which can be the 
case only if the rotation curve rises very quickly in the 
central beam and stays flat further out, as seen the other 
panels of Fig 3. The gas with broad velocity width at the 
nucleus was included when we fit the kinematics and we 
obtained a good fit with a simple flat rotation curve ex- 
tending into the central ring; we did not need to introduce 
any departure from the simple flat-rotation-curve model. 



is 1.0 x 10 6 M Q . The dynamical mass enclosed within the 
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Fig. 2. Top: ATCA first-moment image of H91a + 
H92a emission showing rotation. Beamsize is 1.4" x 1.2". 
Middle: H91cn + H92a rotation curve along the major 
axis of the line-emitting region (PA = 47°). Resolution is 
1.3" x 36kms _1 . Bottom: Rotation curve from the best-fit 
dynamical model described in Sec. 5. 



A. L. Roy et al.: Recombination Lines from NGC 4945 



900 
800 
700 
600 
500 
400 

188 

800 
700 
600 
E 500 
- 400 

f 188 

3> 800 

■| 700 

I 600 
o 

.9 500 
£ 400 



800 
700 
600 
500 
400 
300 



-| 1 1 1 1 1 

. 242 



1 1 1 1 1 1 

244, 



-| 1 1 r~ 



246 



240 



t ■ 



n 1 1 — — i 1 r— 

248 



_j i i i i_ 



i i i i i i_ 



250 



252 



254 



i i i i i i 



256 



V. 



258 



"J 



260 



',262 



264 



266 



268 



.270 



6- 



20 15 10 05 00 15 00 

Minor axis offset / arcsec 

Peak brightness = 3.14mJy/beam 

Contours at 0.8 m Jy/beam * (-3, -2, -1 , 1 , 2, 3, 4, 5, 6, 7 



Fig. 3. H91a + H92a rotation curves along the minor axis direction of the line-emitting region (PA = —43°). Each 
panel is separated by 0.8" along the major axis, with the slice shown in panel 252 passing through the peak of emission 
at the nucleus. The minor axis offset values are relative to an arbitrary offset. Resolution is 1.3" x 36kms _1 . 
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7. Conclusions 

We have discovered H91a + H92a lines in emission in 
NGC 4945 with flux density of 17.8 mJy using the ATCA, 
making NGC 4945 the brightest known extragalactic RRL 
source. The detected line strength infers an ionized gas 
mass of (2 to 6) x 10 5 M and a corresponding star for- 
mation rate of 2 M Q yr _1 to 8 M Q yr _1 depending on the 
model conditions. The star formation rate estimated from 
the RRL detection in NGC 4945 agrees well with rates es- 
timated from radio and FIR luminosities using previously- 
calibrated relations. 

The rotation curve was found to be flat into the central 

1" With Usystomic = 581 km S _1 , ^rotation = 120 km s" 1 . 

We found no need to invoke a bar or radial motion and no 
indication of a kinematically decoupled component. 

Future observations at high frequencies where RRLs 
are stronger and resolution is higher will provide measure- 
ments of multiple transitions to provide better constraints 
on the gas conditions. 

Since RRLs occur over a wide range of wavelengths, a 
given array can provide higher resolution for studying the 
dynamics of ionized gas than can be achieved for studies 
of neutral hydrogen. 

Followup with the Atacama Large 
Millimeter /submillimeter Array of RRLs in NGC 4945 
and others should be rewarding, giving good SNR 
and high spatial resolution to sec finer details in the 
kinematics and of the possible circumnuclear torus. 
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